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^ ■ Abstract 

Q^ , Matrix model approach to multicolor induced QCD based on the 

quenched momentum prescription is presented. It is shown that this 
model exhibits the reduction of spatial degrees of freedom: the parti- 
i tion function is determined by the solution of one dimensional quan- 

^ I tum mechanical problem while the D-dimensional scalar field corre- 

' lators coinside with the same type correlators in the two-dimensional 

induced QCD. 

> 

o 

^ I 1. Investigation of the large number of colors Nc limit of QCD permits a 

deeper insight into the nature of the strong interaction. Analytical sum- 
^ ■ mation of planar diagrams leads to a kind of masterfield equation. The 

. considerably reduction of actual degrees of freedom occurring when Nc ^ oo 

underlies this equation. The large number of " angular- type" variables (uni- 
Qh! tary transformations) is effectively "integrated out" in matrix models. The 

remaining variables form a meanfield described in terms of the masterfield 
equation. 

The case of induced QCD will be considered here. This model was intro- 
duced by Kazakov and Migdal four years ago 0]. Here the induced QCD is 
5^ I revised but instead of the lattice version of QCD in Kazakov-Migdal model 

(KMM) the prescription proposed by Gross and Kitazawa for the planar 
Feynman diagrams is adopted. As in the KMM the absence of the gluon 
kinetic energy enables to integrate over the unitary matrix in closed form. 
However contrary to the KMM case the explicit expression for the "angular" 
integrals are not used in the present approach since it entirely avoids the 
introduction of mean field as well as the solution of the masterfield equation. 
The main goal of it is to establish the connection between 4D induced QCD 
and the same lower dimensional theory. 

2. Since the Gross and Kitazawa prescription plays the crucial role in the 
following treatment it is instructive to begin with brief reminding of it. 

The quantum scalar field in adjoint representation ip{x) is replaced by the 
Nc X Nc matrix if which does not depend on the variables. The field 
derivative is replaced by the commutator 

d^(p --^ i[Pf,, If], 



X 
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where are diagonal Nc x Nc matrixes: 

{p,f = (1) 

The functional (euclidean) integral 

Z = j\{Dv{x)e^v[j d''xTT[--\{d,^Y-\M'v'-V{^))^ (2) 
turns into the matrix integral 

'2tx\° r_ ,2 1 



/rfy,exp(^) {Tr(-[P„y,f--MV^-nV'))} (3) 



where 



dip = Y\.df' 



When Nc oo the diagrams generated by the perturbative expansion of the 
second integral coincide with the planar Feynman diagrams for the first one, 
A being the ultraviolet cutoff. Indeed the matrix model propagator is 

Drawing the Feynman propogator by the doubleline one can assign to each 
line its own momentum pi or p2 so that the total momentum flowing through 
the propagator is k = pi —p2- After this substitution the loop integrals turn 
into the integrals over the momenta pi 2- For planar graphs the momentum 
loops are identical to the color ones inside which the indexes i,j in (|^) circu- 
late. If the D-dimensional hypercube in the momentum space is divided into 
Nc equal cells with the volume /Nc and the vectors in (^) are choosen to 
lay inside these cells the sums over the color indexes in the matrix diagrams 
for the model (|^) yield the integral sums for the momentum loop integrals in 
the planar Feynman diagrams for the functional integral @. When Nc ^ oo 
and A is fixed the sums turn into the integrals. After this the limit A —>■ oo 
can be taken. The hermitean Nc x A^^ matrix (|l]) plays the role of the 
quenched momentum and the components are the values it takes. [| 

The unitary matrixes e*^^ can be treated as a finite dimension approxi- 
mation to the space shift operator. One can define 

V9(x) = e'^'^^e-'^^ 



^It is really not necessary to identify the number of the terms in the integral momentum 
sum with Nc- Indeed, while a planar Feynman diagram contribution is N'^ where 
does not depend on iVc, the rank N matrix model gives for a planar graph Gpianar = N'^ G. 
The quenched momentum prescription ensures G —> G^ for TV — s- cx) and it is this property 
that makes the matrix model to be equivalent to the field theory. Thus only the planarity 
is important here. 
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Then the correlator 

{Trip{xi) ■ ■ ■ ip{xn) ) 

calculated in the matrix model (0) in the leading N^. order will be the same 
as in the scalar field theory (|^). 

The incorporation of gauge field in this prescription is straightforward: 



c 



Here is a Ai'c x Nc hermitean matrix, is a covariant derivative and 
the commutator in the last formula is understood in the matrix sense. The 
(euclidean) matrix action can be written as 

S = TrGl, + 



where Sgi is an appropriate gauge-fixing term. The matrix transformation 
-^A. ^ V~'A^V + -^V-^[P.,V] 

where is a unitary matrix is equivalent through the relations = 
e*^^y4^e~*^^, V{x) = ^^p^Yq-'^p^ to a local gauge transformation in the field 
theory. 

Very important for the following is the additional constraint imposed on 
the measure of the matrix integral, namely, the integration is carried out 
only over the matrixes A^ for which the covariant derivative has the 
same eigenvalues as the matrix P^. An equivalent form of this restriction is 
to rewrite the functional integral as 

\{dA^c^e-^ 



where 

c.r. JdV,6i-^D,-V^-'PM (5) 

and dVfj_ denotes the invariant measure on the SU (Nc) group. If it were not 
for the constraint @ the matrixes would be completely excluded from 
the integral by shifting matrix variables -^^A^^ — > + -^A^ = \D^. The 
constraints have a rather obvious meaning: any component A^ of the 
gluon field can be set to zero by the gauge transformation special for each 
^ (only for the pure gauge field there is a matrix V common to all /i). 
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3. The induced QCD lagrangian describes the theory containing the gluons 
and the adjoint scalars 

£ = \Tr{D,^f - ^TrM'cp' (6) 

but without the pure gluon term. The integration over the scalar field yeilds 
the functional determinant 

Det{M'~D'r'- = exp|-^/iV,logA j d'xTvGl, + Oil/M)^. 

Here 0(1/M) denotes the terms finite when the ultraviolet cutoff A oo. 
They are suppressed by the powers of the scalar mass M which is assumed 
to be very large (although M <^ A). The theory seems to be similar in 
the limit M oo(M <^ A) to the gluodynamic, (7^ ~ 1/ log A/M being the 
coupling constant. One can consider a more general case by adding to the 
lagrangian (||) the scalar field interaction V{ip) = J^'^nf"'- 

According to the quenched momentum prescription the matrix model in- 
tegral for this theory in the hamiltonian gauge = reads 

Y[dA^ c^,6{Ao)dip 



X exp(^)"rr|lB^.>^f-^^^ 



Resolving 5-functions in the constraints gives 



Z = JlldV,d^6(^^ (Po - V,-'PoVo)^ 



Writing here the scalar field matrix as = f/ ^(^dU where f/ is a unitary 
matrix and ipd is a diagonal one the integral takes a form 



Z = [HdV^dU difd A\^d) 5 — (Po - V^'PoVo 

V{ipd) = Y.iV{ipi) is the sum over the eigenvalues, A{ipd) is the Van der 
Monde determinant. Changing the integration variables 

V,U-' = V„ ix=l,...,D-l, (7) 
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allows to separate out the gauge-fixing term: 

Z = J dVo6 (Po - K)"'^oK))) ■ / dUYl dV^d^aA\^d) 

Dropping the tildes and redefining Uq = Vq the partition function can be 
rewritten as (a ~ l/A^(Po) is the normalization constant) 




The integration over the unitary matrixes in decays into the product of D 
equal integrals (they are equal since the matrixes P^ are unitary equivalent). 
Consider the integral 

Z^ = JdV exp (^^y [^Tr[P, , Vy^.V-']'] (9) 

for the fixed matrix ipd and index fi. An explicit calculation of this expression 
is a rather nontrivial problem although an essential progress has been made 
recently in studying the similar type integrals 0. Instead calculating it the 
much more weak "scaling" property of this integral will be enough for the 
following. 

The main contribution to the integral over the matrixes if in (^) comes 
from the domain where (pij ~ 1 and where the typical scale of the eigenvalues 

~ ^y^- It is just the region that gives the leading contribution of the 
order A*"^ to the free energy while the measure of the rest integration domain 
tends to zero in the Nc ^ oo limit. In this limit the eigenvalues distribution 
is described by the smooth density function p(A). 

Since Tr[P^, v?]^ ~ one can expect 

= exp{N^^F + 0{N,)} (10) 

where the coefficient F depends on the eigenvalues of the matrix (f and for the 
large can be treated as a functional F = F[p]. Eq.(|lOD is a consequence 
of the classical nature of the large Nc limit which manifests itself in the 
factorization of the correlators of colorless {U{Nc) invariant) operators 0|. 
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The large Nc behavior of the unitary matrix integral (|^) looks like as if 
it is dominated by a saddle point. There is a representation which makes 
the dependence in the integrand to be explicit 0. In the large Nc 
limit SU {Nc) algebra is equivalent to the infinite dimensional Lie algebra 
of area preserving diffeomorphisms of the sphere SDif f{S'^). The matrix 
A from SU{Nc) algebra transforms in this representation into the function 
A{xi,X2,X3) on the unit sphere x\ + + x\ = 1, the commutator being 
replaced by Poisson braket 

hm ^[A,B] = {A,B} 

which are defined as 

{A, B) = XiEiki T— T—- 
dxk dxi 

The integral (|^) turns into the functional (infinite order) integral over the 
group of area preserving diffeomorphisms 



Z^ = j Wexp \^Nl J dQ{p, V{v)} 



(11) 



where V{ip) denotes the action of the diffeomorphism on the function cp and 
/ dQ is the integral over the unit sphere. The functions (p and p correspond 
to the matrixes (p and in the integral (|^). Their particular structure as 
well as the result of the action of V on ip is not important here. Only one 
thing is significant for the following, namely, comparing the expression (pT|) 
with (0) one can conclude that 

JdVexpi^-^a (^^y Tr[P,V-'vDVf^ = exp {a iV^F + 0(iV,)} 

since it is nothing more than rescaling of Nc by the factor ^/a. Thus the 
power of the integral can be replaced within the leading order accuracy 
by the single integral 

{Z^f = JdVexpS^-^D (^^y Tr[P,V-'^dV]'y (12) 

The " scaling property" ( |T^ ) will be the central point for the following treat- 
ment. All other results are more or less trivial consequences of it. 
The relation (|T2p applied to the integral (|) gives 

Z = a j dipd^^{Vd)dV 

X exp {^-^y Tri^-D \P, , V^dV-'f - - V{^d)\ ■ 
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After combining the factors in the integration measure the partition function 
M) takes the final form 



Z = ajd^exp (^^y Tr j^D [P^ - ^MV^ - ^(v^.)} • (13) 

Here the matrix index is fixed by the value = 1, the result being clearly 
independent on the particular choice of the space direction. 

The expression (|T3|) is a Gauss integral for V{ip) = 0. In the leading 
order 

logZ = const + DNl log A - ]^Y.^og[D{p\ - plf + M^] = 

ik 

= const + DNl \ogK- -NIVd(—] / / \og{Dp^ + M^) 

where Vr, is a total space volume^. Collecting the terms independent of M 
in the factor Sq one gets for the vacuum energy density 

4. When the scalar interaction V{}p) 7^ the evaluation of the integral 
is equivalent to solving the quantum mechanical problem. Indeed, as is seen 
from inverting the quenched momentum prescription, it coincides in the large 
limit with the partition function 

Zd = y'n^</'(^) expTry'rf^x|-^D(9i(/.)2-lMV'-V^(<^)} (15) 

Since there is only one actual degree of freedom contributing to (|l^) 



A 



where E is the ground state energy for the one-dimensional system described 
by the integral 

Zx = JllVQ{t)expN,TrJdtl^-^Di^Q^-^l2M^Q^-V{Q)'j (16) 

in which fi = A/2n and V{q) = J2n>3 f^^~^^^"'^'^~^^Vnq"'- This relation immedi- 
ately follows from the lattice version of the theory (plsf ) which is equivalent to 
the ultraviolet regularization with the cutoff A, being the lattice spacing. 

^The value Vtot = {2tt)'^ Nd is to be taken as a space volume in the quenched model 
prescription to recover the planar graphs. Taking into account the proper order of limits 
(first Nc — > 00 and then A — > 00) one can put Vtot = NcVd 
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The system (|T6|) is solved by reduction to the free fermions moving in the 
external potential V{q) 0. The energy is E = Yjf^iEi, where £i are the 
energies of the lowest Nc occupied states: 



+ - /iM^ + V{x) 



ipi{x) = Siipiix). 



2fiD dxl 2 

In the quasiclassical approximation which validity is justified by large A^, 

dpdq 



E = Nt 



277 



H{p,q)e{S-H{p,q)) 



(17) 



where 9 is the step function, 

1 

2D ' 2" 



n>3 



and Fermi-level S is determined by the equation 

dpdq 



J Zn 



Note that E = for M = and V{q) = 0. 

5. Consider now the two-point correlator of the scalar fields: 



= ^ Jl[dA,c,dvTr[^{xMO)]HAo) 

1-1=0 

X exp(^)%r|lx:[^.,^f-^^V'-n^) 
The variables changing (0) brings the integral to the form 



(18) 



1 



D~l 



dVo n dV^dUd^dA\ipd)S 



Nr. 



Po - V.-'PoVc 



(19) 



X exp 



27r\^ „ rl 



A 

D-l 



11=1 J 

In contrast to the partition function case the pre-exponential factor prevents 
to factorize the gauge- fixing term with the integral over Vq. Nevertheless the 
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"dimensional scaling property" (|12|) can be used again to reduce the product 

■D 



of the independent integrals ^ to the single one: 



K^^\x) = \j dVodVidUdipdA^iif,)5 ^ (Po - V.^'PoVq 



X exp^- 



2 

Now the steps leading from (|T^) to ( ]T^ can be repeated in the inverse order 
which yeilds 



1-1=0 

X exp (^^y Tr {^[Do, + \{D - 1) [D,, 



After rescaling 



Po VD^P„ (20) 

" - 7F3T^*'V (21) 

the integral takes a more symmetric form | 



X exp ( - j Tr <j - V? - 7^7^ - V{^) 



where 

= ^(D-l)-t/i(f-i)("-2)t;„y," (22) 

n>3 

or 

ir(^)(Xo,Xi,0,...,0) = ^^-2_1_K(2)(v/:D^Xi,X2). 



■^The first line here is the quenched momentum analog of xq xq^/D — 1 replacement. 
Although the integration region is no more a hypercube in the momentum space it is not 
essential for a ultraviolet convergent 2D theory. 
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Here K^'^^ is the correlator in two-dimensional theory. 

Thus the two-point scalar correlator in the D- dimensional induced QCD 
coincides with the same (up to rescaling) correlator calculated in the two- 
dimensional induced QCD. This result clearly holds for any ra-point scalar 
correlator provided all the points lay in the same two-dimensional plane. 
6. Unfortunately the above methods based on the "scaling property" (0) 
are insufficient for gluon correlators. The tensor structure of the vector field 
correlators in D-dimensional space is much more rich than for D = 2 and 
they can not be reduced to a two-dimensional theory. However there is a 
special kinematic to which a little modification of the property (0) can be 
applied. 

Consider again a two-point correlator 

d^^'^x) = {Trn^A^{x)n,A,{0)) (23) 

where n is an arbitrary vector of the form = {0,ni, . . . ,nu_i) (recall 
the gauge is Aq = 0) normalized as = (D — 1). One can always choose 
the coordinate axes in such a way that = (0, 1, ... , 1). The vector x in 
will be assumed to be of the form x = {xq, Xi, 0, . . . , 0) in this basis, so 



ixn)'^ = Xi = x\ 



It is convenient to deal with the generating functional which in the quenched 
momentum prescription is 



j\idA^c^5{A^)d^ (24) 



D-l 

n 

/i=0 

X exp(^)%r|lx:[^.^V^f-^^V'-n^) 

^^ n ) 

The correlator is clearly given by differentiating ( p4D with respect to the 
matrixes Ji 2'- 
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Jl,2=0 

After variables shifting and resolving the constraints Zj takes the form 



exp|^rr(J2 + Ji) 

exp (f)''r-{^[^o, -\m^^- -V{^) 
jlidV, exp (^)'' [Tr\[V^'P,V,^^f 



(25) 
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The internal integral here is again the product of {D — 1) equal "angular" 
integrals. Although they differ from by additional terms in the ex- 

ponent these terms are of the order A*"^ {Jik ~ 1) therefore the same type 
saddle point behavior is natural to be assumed for these integrals too. The 
continuous representation like ( [TI| ) can be also written for each of them as 

27r\^ r _ ri 



where V{p) denotes as before the action of the area preserving diffeomorphism 
V on the continuous image (function on the unit sphere) of the matrix P^, 
J2{x) and Ji are the continuous images of the matrixes e~*^^J2e*^^ and Ji 
respectively. This formula enables to replace the product of the integrals in 
( |25D by the single integral for the momentum Pi. Introducing then the field 
Ai through the relation 

v-'p,v = -Di = Pi + ^/if-Mi 



and making rescaling (pO]), (pT]) the generating functional takes a final form 

Zj = Nexply^Tr{J, + J,)J2P^ - {D-1)pA 
[9 fi ) 

J n dA, c, 6{Ao) exp Tr | i Y.[D,, v] 



1 M2 2 

-V{^) (Z}-l)/xf+V2Ai(x) - (D-l)^T+ViAi} 

where ^(v?) is given by (^) and the factor does not depend on Ji^2- It 
gives for the correlator 

d(^)(xo, xi, 0, . . . , 0) = (D - 1) (TrAi(v/D^a;o, xi) Ai(0))2 (26) 

where the right hand side means the correlator in the two-dimensional in- 
duced QCD. 

Note that a two-point gluon correlator in the dimension D 

dlS\x) = {TrA,{x)AM) 
is determined by two scalar functions 

dj,u\x) = xf^x^ d^^\xo,xl) + 6^^di^\xo,xl). 
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The formula ( p6D implies one relation between them: 

= {D- 1) V' [di'\VD^Xo,xl) + d!i'\^/D^xo,xl)]. 

There is another way to derive this relation. Instead of ( p3D one can start 
from the correlator 

dS\x) = Y.{TrA,{x)A,m = ^{Tr[P,P^ - D,{x)D,m) 
fi 9 

and introduce the generating functional 

= {exp{aTrD^ix)D^{0)}), 

(7=0 

The integral for Z^ allows for the dimensional reduction like (^4|) since the 
"angular" part of it is again the product of {D — 1) independent equal inte- 
grals. 

One can combine both these trick to derive similar relations for higher 
correlators with the gluon fields. However there is a lot of tensor structures 
which are not involved in such relations at all. For instance, the triple gluon 
vertex does not contribute to any of them because of antisymmetry over the 
index permutations. 

7. Essential reduction of the space degrees of freedom occurring in the in- 
duced QCD is the main result of this paper. All D-dimensional induced 
QCD theories turn out to be related through the equality ([T^). The reason 
for the "scaling property" (|12D lies, probably, in the fact that coordinate 
and momentum operators can be approximated by large N^. matrixes. Al- 
though there are no finite order representation for the Heizenberg algebra 
one can construct the matrixes and from SU{Nc) algebra for which 
[X^,P[/] = iSfj^i, + 0{1/Nc) It explains why space disappears in the in- 
duced QCD: the space-time transformation can be mapped for any dimension 
D into the same SU (Nc) group. The quenched momentum prescription is a 
particular choice of such a map. One should stress the point however that 
the dimensional scaling is valid only for the gauge theories because only for 
them the "angular" (over unitary matrixes) integrals are carried out sepa- 
rately for each field component. It would be impossible for a scalar theory 
where the unitary rotation is common for all space directions. It is the ad- 
ditional unitary integrals in the QCD that "absorb" the space degrees of 
freedom. 

Thus the induced QCD results into the low dimensional theory - one di- 
mensional quantum mechanics for the partition function case and two dimen- 
sional theory for the scalar correlators. The latter fact is of special interest 
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because the two dimensional QCD models with a scalar matter in adjoint 
representation have received a recent attention due to a similarity they have 
with 4D gauge theory [|, The scalar degrees of freedom "compensate" 
the absence of transverse gluons in 2D models. From the other hand it is 
a commonplace that the linear potential resulting from the Coulomb inter- 
action in 2D theories has no relation to 4D confinement. The dimensional 
reduction hints that such a relation might exist. Indeed, the interaction po- 
tential between the large mass matter particles can be extracted from the 
amplitude of their elastic 2 — >^ 2 scattering. It is the amplitude that allows 
for the two-dimensional reduction in the induced QCD. The Coulomb force 
originates from term in the effective lagrangian which appears after inte- 
gration over the scalar field. However there is not the logarithmic divergency 
when D = 2, and the first term does not dominate in the loop expansion 
in which the higher terms are related to the presence of the non-Coulomb 
degrees of freedom. 

The natural question here is in what extent these results are valid for 
the real, not induced, QCD. Unfortunately the naive picture supposing all 
the terms except the first one in the loop expansion to be suppressed by the 
large scalar mass M is inconsistent. The scalar loops result into the point-like 
gluon interaction only if the gluons' momenta are much smaller than the M. 
Therefore the induced QCD has to imply two typical scales to be equivalent 
to the gluodynamics - the ultraviolet cutoff A for the scalar particles and 
their mass M -C A for the gluons. However the gluons' momenta circulating 
in the internal loops are restricted really by the A rather than the M value. 
That is why the expressions (|I^) , (|T7p for the vacuum energy density can not 
be immediately applied to the QCD since the hard gluon loops contribution 
is not separated from them. 

A possible way to connect the induced QCD with the real gluodynamics 
is to adopt a slightly different approach. One can start from the scalar QCD 
treating the scalars on an equal footing with the gluons so the lagrangian will 
encounter the scalars selfinteraction hip'^. The constant h as well as the gluon 
coupling g are the bare constants which have to be taken to be functions of the 
regularization parameter A. They values are determined by regularization 
procedure: the gluon and scalar dressed vertexes and the residues of the 
propagators poles are fixed at some external momentum scale /i referred to 
as a normalization point. The position of the scalar propagator pole is fixed 
too as a physical scalar mass. If the bare constants are adjusted to keep 
all the fixed at the point values when the ultraviolet cutoff A cxd, then 
according to the renormalization theory the Green functions will be finite. 
It is the dressed vertexes values at the normalization point rather than the 
bare constants that play the physical charges role in the theory. 

The bare gluon coupling is zero in the induced QCD, therefore the total 
number of the bare parameters is insufficient to satisfy all the normalization 
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conditions. However one can lift one of them, say, not to fix the position of 
the scalar propagator poles using the bare scalar mass as a parameter to keep 
the gluon vertex value. Suppose the renormalization scheme can be carried 
out for all the vertexes and residues while the physical scalar mass turns out 
to be of the order of A. In this case the renormalized theory will be surely 
equivalent to the QCD. Even if the more complicated situation occurs it will 
be of interest to study what is the limiting theory. The obtained here results 
enable to investigate this problem in the framework of the two-dimensional 
scalar QCD. 

There are several difficulties however to proceed in this manner. It is 
unclear if the triple gluon vertex allows to be reduced to the two-dimensional 
theory like the propagators do. The possible way to circumvent this obstacle 
is to determine the gluon physical coupling through the interaction with a 
large massive matter field, that is through the Coulomb potential. An other 
probably more serious difficulty is the absence of an exact solution of the 
massive two-dimensional QCD although there are numerical investigations 
of it i. 

Note in the conclusion that the induced QCD could be interesting in itself 
as a toy model even without the direct equivalence with the gluodynamics. 
Indeed it provides an example of the four- dimensional theory asymptotically 
free at the small distances (in some sense the asymptotic freedom in the 
induced QCD is even more strong than in the gluodynamics |jlO|) which 
exhibits a confinement-like behavior at the large ones. 
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